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The helical absolute equilibrium of a compressible adiabatic flow presents not only the polar¬ 
ization between the two purely helical modes of opposite chiralities but also that between the 
vortical and acoustic modes, deviating from the equipartition predicted by Kraichnan, R. H. 
[1955 The Journal of the Acoustical Society of America 27, 438-441.]. Due to the existence of 
the acoustic mode, even if all Fourier modes of one chiral sector in the sharpened Helmholtz 
decomposition [MoSES, H. E. 1971 SIAM (Soc. Ind. Appl. Math.) J. Appl. Math. 21, 114-130] 
are thoroughly truncated, leaving the system with positive definite helicity and energy, negative 
temperature and the corresponding large-scale concentration of vortical modes are not allowed, 
unlike the incompressible case. 
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1. Introduction 

“Local anisotropy’ or “polarization”, referring to the dependence on the directions in the local 
coordinate frame of the w ave vector space, is one of the important characteristics of turbulence 
( Cambon & Jacqui'nlll989i and references therein). Such polarization can be due to either linear 
(damping or wave) or nonlinear mechanisms (helicity, s ay), or both, as appears in ro t ating flows 
and qu asi-static magnetohydrodynamics, among others dCambon & Godeferdlll993l) . Kraichnan 
( I 955 I ‘K55’ hereafter) studied the absolute equilibrium (AE) of an adiabatic compressible flow 
model and showed with small-excitation approximation the complete equipartition of energy 
among all the vortical , dila t ationa l and density mode s [similar result for the incompressible case 
had been obtained by lLeg(ll952h and iHopf I (Il952h . the latter of which applied functional cal¬ 
culus to formally deriye it, but without explicitly introducing the ne cessary Galerkin t runca tion], 
and non-equipartition was anticipated for the strong-excitation case. Krstuloyic et al.l ( 2009h pro¬ 
posed a new algorithm to generate the irrotational, thus nonhelical, compressible AE, and they 
indeed obseryed nearly Gaussian and equipartition at lower temperatures but non-Gaussian statis¬ 
tics and non-equipartitions of total compressible energies, together with “non-physical” negatiye 
fluid density for higher temperatures. K55 also discussed th e effects of po larized dissipation. 
Note that helicity is (formally) conseryed in barotropic flows ( Moffatdl 19691). thus, acco rding to 
Kraichna^ ( 1973h . incompleteness of K55’s analysis. Also note that KraichnanI ( 19731) refuted 
the speculation of the possibility of inyerse cascade in three-dimensional (3D) helical turbulence 
from the superficial analogy with 2D flows bearing also two ideal inyariants, which may lead 
to the belief that two positiye-definite inyariants must indicate inyerse cascade: We will update 
K55 by additional respect of helicity and show that the acoustic mode in the compressible flow 
excludes such a possibility. We should notice that such multiple-constraint analysis based on 
time-rey ersible system cou ld haye releyant indications for nonequilibrium and irreyersible en- 
sembles ( [Gallayotti 2014 ) of other types of time-reyersible modified Nayier-Stokes equations 
dShe & Jackson Ill993h as we will come back in the end. 
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Let u, p, 7] and c represent the flow velocity, density, viscosity and velocity of sound in a 
compressible flow with the pressure p given by the adiabatic relation p = c^p. Since now, p is 
independent of the temperature (otherwise the temperature should be fixed, i.e., iso-temperature), 
the mass and momentum conservation laws are decoupled from t he energy equation, formi ng an 
autonomous system where the entropy modes are not concerned ( Sagaut & Cambon 2008h . K55 
introduced the transformation p = poexp{C} (‘Kraichnan transformation’ hereafter) with pq 
being the equilibrium density, which, in the forceless case, leads to the equations of motion 


C “f “f — 0 

“f C ^^A ^^A(T,cr — 0: 


( 1 . 1 ) 

( 1 . 2 ) 


where 9ap = Ua^p + up^a — ^^‘pUa,(T, (•),7 = d{»)/dx'^ and (•) = dt{»). Working in a cyclic 
box of dimension 27r with V = [0, 27r)^ and applying Fourier representation for all the dynamical 
variables v(r) v(k), say, u(r) = u(k) exp{ik ■ r} with = —1, K55 constructed a 
phase space by the real and imaginary parts of vs and showed that the flow in such a phase 
space is incompressible, the Liouville theorem, in the inviscid case. K55 further took the ‘small- 
excitation’ approximation for the ideal invariant of the total, kinetic plus potential, mean energy 
per unit mass: 


^=2 

k 


(1.3) 


which, according to Eqs. (Ol and (11.2b with the viscous terms removed, obeys in r-space 
E + (uaC^a) + {{(ua),a) — 0 with 6 = UaUa + and (•) = ( 2 ^ Iv The peri¬ 

odic boundary condition kills {{Cua),a), and E is not rigorously an invariant but Kraichnan’s 
approximation of the total energy. Galerkin truncation, say, imposing all modes with |fe| greater 
than some cut-off value K to be zero, is then performed, which does not c hange the approxi - 
mated conservation of energy due to its being quadratic and diagonal in k (lKraichnanlll973h . 
Reasoning with the iF-theorem, K55 predicted that an ensemble of systems will tend towards 
an equilibrium distribution uniform over the surfaces of constant energy, thus an equipartition, 
in particular between the shear/vortical and acoustic modes: Such K55 fields are isotropic, both 
globally/‘directionally’ (independent on the direction of fc) and ‘locally’ (independent on the di¬ 
rection s in the local coordinate frame), in the terminology of Cambon & Jacauinl(ll98^ . lMoffattl 
( 19691) has shown that, with the vorticity defined by u;(r) = V x u(r) = ‘^{k) exp{ifc ■ r}, 

the helicity 


H= - ^ M(fe) • d>(-fe) 


(1.4) 


is also fo rmally an invarian t in barotropic flows with p = p{p) and, just as in the incompress¬ 
ible case ( Kraichn^ 1973 ). should be r espected. It is also desirable to extend the polarization 
analysis to cases with strong excitations ( Shivamoggill 9971: iKrstulo vie et al.ll2009ll . with ad hoc 
assumption or uncontrollable approximation, but we will only focus on the polarization issue, 
using still the small-excitation approximation for preciseness and analytical tractability. Our pur¬ 
pose is to deliver a crystal clear message of the so-called ‘helicity effect’ in such a situation, 
instead of some heroic theory or the final solution of the whole problem (with general equation 
of state and strong excitations) which appears extremely complex, as we will further remark later: 
indeed, since the absolute equilibrium is nevertheless far from the turbulent state (except for the 
large scales in the inverse cascade case) and in general is for qualitative theoretical insights, it 
would hardly be more illuminating for our purpose to further write down or compute the very 
complicated strong-excitation/high-temperature absolute equilibria. 
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2. Compressible helical absolute-equilibrium polarization and non-equipartition: 
differences with the incompressible case 


MosesI (119711) has ‘sharpened’ the Helmholtz theorem by further decomposing the transverse 


(‘solenoidal’/'vortical’) velocity field into a left-handed and right-handed chiral modes of sign- 
definite helixity (helicity intensity) 


u(k) = u+{k)h^{k) + u_{k)h_{k) + u^{k)k/k 


( 2 . 1 ) 


with ik X hs = skhg and s = ± (denoting opposite — right- v.s. left-handed — screwing 
directions, or chiralities, around k), and that 

i? = - ^ |w+P + (2.2) 

k 

for small excitation as in K55, and 

( 2 . 3 ) 

k 


The equations for incompressible flows (with out the longitudinal/dilation al c omponen t s) cor- 
responding to Eqs. (12.112.212.31 l were used by Cambon & JacauinI (Il989l) and IWaleffel d 19921) 
to study rotating turbulence and triad instabilities etc. Denoting u±{k) = u+{k)h^{k) -f 
U-{k)h^ (fc), we can rewrite Eqs. (II.Il l and (11.21) without viscosity in the following way 


'm= E [l{p)p-u±{q) + Q{p)p-u\{q)\ + ku\{k) 

p+q=k 

ft|(fc) = N^ik) + c^kCik) 
ii±{k) = N^{k), 


(2.4) 

(2.5) 

( 2 . 6 ) 


with N = u ■ Vu = (it| -t- u±) • V(t6| + u±), showing the complicated interactions [with, 
say, iV|(fc) = (mj_ • Vm_l)| + {u\ ■ Vm|)| + (mj_ • Vm|)| -f {u\ ■ Vm_l)|], of the acoustic and 

vorticity modes. Note that the terms C{p)p ' u±{q) and {u± ■ Vu±)^ as in the incompressible 
case conserve, for each interacting triad, and (u^) respectively. Small-excitation approxi¬ 
mation for t he total energ y in terms of C, and u is indeed necessary to proceed with the analytical 
calculation. MosesI ( 197 ih solved the corresponding linearized problem, and K55’s treatment can 
be viewed as quasi-linear [with = p'{po) being the constant of the background field The 
inviscid and forceless Galerkin truncated dynamics for absolute equilibrium analysis may also 
be viewed with the balance of the nonvanishing internal viscous force and the external stirring 
force f{t). It is trivial to take the external force be the internal viscous force with constant vis¬ 
cosity q, but it may be more interesting to modify the viscous term to balance the independent 
external forcing; One particular possibility is, say, replacing the q neglected in Eq. (I2.51 l for the 
longitudinal component with a dynamical eddy viscosity p| (k, t), i.e.. 


p-^P|(fc) = 3/|(fc) •«! (fc)/[4fc2|M|(fc)|" 


(2.7) 


and similarly q± for the left and right chiral sectors of Eq. (I2.61 l f or the transverse compon ents, 
mutatis mutandis. Eollowing the nowadays a routine procedure (IZhu. Yang & Zhull2014l) . we 
immediately obtain with the constant of motion C = aE + flH from the Gibbs distribution 

























J.-Z. Zhu 


exp{—C} the absolute equilibrium modal spectra 

U+{k)^{\u+\^) = 

U-{k)^{\u_\^) = 
Z(fe)4(c2|Cp) = 
C/I(fe)4(|«||2) = 


1 


a + /3fc ’ 
1 

a-13k' 
1 

a 

1 

a 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


We see that the spectra are directionally isotropic; but, when (3 or H is not zero energies are 
polarized, except for the equipartition between the dilatational and density/pressure modes; and, 
the vortical-mode energy is larger than the ‘acoustic’ one 

U-^ = U++ U->U~ = Z+ U^, (2.12) 

which indicates that helical states tend to have higher vorticity-mode energy and reduce ‘noise’; 
see more discussions in the next section. [It is interesting, though probably useless in practice, 
to note that locally in the plane perpendicular to k, the statistics are still isotropic, because the 
polarization between U'^ and U~ is about the screwing directions around k, having no pre¬ 
ferred direction perpendicular to it.] A physical question is whether and how such polarization 
and equipartition of the thermalized state persists in the nonequilibrium turbulence, in particular 
whether the acoustic-mode energy level stays in between the vortical ones of opposite chiralities. 

Note that Eqs. (12.8112.91 1 are respectively for ‘pur ely helical absolute equilibrium’ of the ‘chi- 
roids’ Ushs{k)e*‘''^ + c.c. (IZhu. Yang & Zhull2014l) . Due to the existence of the acoustic modes 
(u| and Q, even if we thoroughly remov e one chirality, say, that denoted by of the vortical 
modes, as in the numeri cal simulations ofiBif erale, Musacchio & Toschi ( 2013h and in the math¬ 
ematical proof used in Biferale & Titi (Eolsi) . realizability of acoustic modes (positiveness of 
C/I and Z) exclu des the possibility of ‘neg ative temperature’, a < 0, unlike the incompressible 
homochiral case ( Zhu, Yang & Zhull2014 ): so, there is generally no indication, from such abso¬ 
lute equilibrium, of large-scale concentration or inverse transfer of U for the compressible case, 
except that the large-scale concentration state with the very specific truncation scheme discussed 
in Sec. 2.2.2 of Zhu, Yang & Zhul ( 2014 ) may stil l be possible. On e thus m ay say that helicity has 
less pronounced effects for compressible flows. Krstulovic et (I2nn9h instead further implic¬ 
itly truncated the rotational modes by using the equations of irrotational flow [otherwise, since 
the vorticity is not invariant after Galerkin truncation of the original Euler equations, it would 
not make much sense to talk about compressible irrotational absolute equilibrium, even if with 
an irrotational initial field, as they did]: In the numerical results without limiting to the small- 
excitation-approximation formulation, nearly equipartition at low temperature and obvious non- 
equipartition and non-Gaussian at high temperature were observed; the formulation they used 
for more general, but still isoentropic and barot ropic, flows a pplies the transformations of We¬ 
ber and Clebsch to find the variational principle ( Mobbsl[l982 ). which involves several auxiliary 
potential functions, making the problem structurally complicated and practically inconvenient. It 
may nevertheless possible to proceed the calculation and computation, though formidable and 
difficult to be il luminating [even for th e very specific irrotational case, the calculations and 
computations of iKrstulovic et all (2009) offer no new physical insight except for the general 
notion of nonequipartition/non-Gaussianity already expected by K55]. Actually, based on the 
above fundamental results, various other approache s can be further app lied for different specific 
problems: Eor example, the deci mation sche mes bv Frisch et al.l (120121) and the variety of chiral 
(sub)ensembles in the fashion of Zhu ( 20141) can also be constructed for particular problems. 

The effect of the Galerkin truncation on a single solution can be well understood in the one- 
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dimensional hyperdiffusion and hyperBurgers model as shown by Bovdl ( 1994I) . who resem¬ 
bles the “infinite shock”, as the asymptotic Galerkin truncation from hyperviscosity of spec¬ 
tral viscosity, and the “Gibbs shock” (the asymptotic approximation to the truncated Fourier 
series of a discontinuous function, the ‘Gibbs phenomenon’). Note that it defies analytical solu¬ 
tion to the Burgers models to intr oduce Galerkin truncation and hyperviscosity or other general 
spectral viscosity ( Tadmoi 1989ll . Indeed, Galerkin truncat ed Burgers has intrinsic stochastic - 
ity which will approach the absolute equilibrium in the end dAbramov, Kovacic & Maid^l2003h . 
and the residue of thermalization leads to the bottleneck phenomena in the hyper viscous sim¬ 


ulatio n , like that of Navier-Stoke s, in between the inertial and dissipation ranges (IFrisch et al. 


2008 ). Sagaut & German^ ( 2005h observed that filtering a shock introduces parasitic contribu¬ 
tions which can seriously spoil the small-sca le turbulence. Neverthel ess, various fingerprints 
of thermalization may be found in turbulence ( Zhu, Yang & Zhu 20l4 and references therein), 
suggesting the potential of helical absolute equilibrium for understanding shock dynamics in 
the possibly helical three-dimensional Burgers’ turbulence (c ^ 0), generalized for rotational u 
(which in the traditional definition is the gradient of a potential), or compressible turbulence with 
pressure. 


3. Effects of polarized thermalization and dissipation in turbulence 

The dissipation rate for the energy of the transverse-velocity (u±), longitudinal-velocity (u|) 
and density (Q modes are different, respectively rjk^, and 0 (neglecting the ‘second’/‘bulk’ 
viscosity). Besides the small-excitation assumption, K55 further assumes the modes “exchange 
energy very rapidly” and concludes that the ‘average’ dissipation rate of ‘turbulent [solenoidal/vorticity- 
mode] energy’, Tjk'^, be larger than that of ‘acoustic energy’, ^pk"^, thus slower decay of turbu¬ 
lence with ‘noise’ than that without ‘noise’. The theoretical assumption concerning time scales 
should better be related to the spacial scales/wavelengths (see below). Note however that the 
smaller ‘average’ dissipation rate of acoustic modes is due to zero dissipation of the density 
modes, and the compressive velocity modes themselves have larger dissipation rate than that 
of the vorticity rn odes: > pk"^. It is interesting to examine the simulation results by 

Wang et al.l (l2013h . though it is not completely clear how much relevant are the adiabatic and 
small-excitation assumptions made in the calculation; In their simulations the total compressive¬ 
mode energy is much larger than the vortical-mode one, which however happens in the largest 
scales and which might be due to the forcing mechanism (they forced on both modes but did 
not tell further details which can also enter to shape the spectra); while in the inertial range 
with well-developed turbulence it is shown that the compressive-mode energy level is lower 
than that of the vortical modes, with apparently steeper scaling law. We now try to understand 
these observations by the polarized thermalization and the polarized dissipation. The time scales 
between the quadratic and linear terms in Eq. (12.41) and (12.5b balanced (‘crit ical balance’ as¬ 
sumption as in magnetohydrodynamic turbulence ( Goldreich & Sridhan 19951) as our tentative 
step for theoretical interpretation), the eddy turnover time and/or the time scale of acoustic-mode 
interaction/conversion can be estimated to be l/(c/c), while the dissipative time scale for the 
compressive mode U| is {^pk^)~^. When the dissipative time scale enters and 




(3.1) 


the internal exchange between the density and compressive modes is faster than dissipation, thus 
K55’s ‘av erage’-dissipation-r ate argument works, with slower decay of the acoustic modes as 
appears in IWang et al.l (l2013h ’s simulation with normal viscosity; while in Wang et al.’s simula¬ 
tion with hyperviscosity, with the exponent n = 8 in the dissipation quickly becomes faster 
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than the internal conversion with 




(3.2) 


as k increases, thus K55’s ‘average’-dissipation-rate argument won’t work, resulting in faster 
decay of the energy spectrum of the compressive mode than the solenoidal mode. If the flow 
is strongly helical, then polarized thermalization effect may enter: in the former case with nor¬ 
mal viscosity, vortical mode should be stronger with the decay slowed down, competing with 
polarized dissipation effect; in the latter case with hyperviscosity and when K55’s argument 
does not work, both the polarized dissipation and the polarized thermalization favor the vortical 
mode, leading to the strengthened de parture between the vortical and dilatational mode. Thus, 
the speculation by Wang et al. ( 2013h that “As Mach number increases, the interactions between 
the solenoidal and compressive modes will become stronger, probably leading to a tendency of 
energy equipartition between the two modes”, suggested by K55 as the authors explicitly re¬ 
ferred to, may not be as favorable for strongly helical compressible turbulence: our Eq. (12.121 1. 
t/-*- > t/"^, in the absolute equilibrium state might be extrapolated to more general cases, which 
could be preliminarily seen by perturbing the last term of Eq. (I2.21 i and simply check the cor¬ 
responding results. Unfortunately, it seems to us that no documentation of compressible helical 
turbulence data has been made available, and relevant in-silico or laboratory experiments are 
ca lled for. __ 

Bertoglio, Bataille & Marion ( 200lh reported in the two-point closure theory (with also the 
extra small-excitation assumption, represented in the low Mach number and the neglection of 
density fluctuations or entropy modes etc., in their “EDQNM” approach) and simulations, forced 
only on the solenoidal modes, asymptotic equipartition between the kinetic energy and the po¬ 
tential energy of the acoustic (‘compressible’) mode (Eig. 17 there), which is consistent with 
the absolute equilibrium equipartition, suggesting that K55’s assumption that the acoustic modes 
‘exchange energy very rapidly’ may work in their model and simulation so that the dissipative po¬ 
larization is not effective. However, the EDQNM closure involves approximations or ad hoc as¬ 
sumptions, especially as pr oposed in the subsequent st udy by Fauchet and Bertoglio [according to 
the detailed discussions by Sagaut & Cambon ( 20081) 1. more careful treatment for the renormal¬ 
ization of the dispersion frequency of acoustic waves may be necessary and can lead to the break¬ 
down of the so-called “strong acoustic equilibrium”, the equipartition between the dilatation and 
density modes at eac h k as appears in our stati stical absolute equilibrium spectra and in the purely 
linear-wave regime ( Sagaut & Cambonll2008h . Absolute equilibrium concerns the ultimate state 
of thermalization whose times scale is not necessarily much smaller than the turbulence eddy- 
turn-over time scale, so it would not be surprising that such precise “strong acoustic equilibrium” 
break down in the turbulent state as a result of delicate balance among various mechanisms with 
different time scales. If we relate the turbulence eddy turn-over time to the de-correlation time 
modele d by the eddy-damping de-correlati on function, then indeed the traditional exponential 
one by iBertoglio . Bataille & Mari onl (1200 ih is much slow er than the Gaussian one by Fauche t 
and Bertoglio ( Sagaut & CambonI 20081) . so the model of Bertoglio, Bataille & Marion ( 2001 ) 
may allow enough thermalization time to have the strong acoustic equilibrium. Note that we have 
related the strong acoustic equilibri um state to the statistical equilibria, besides the classical lin¬ 
ear waves ( Sagaut & Cambonl2008 ). So, our results could provi de insights for practical us e, such 
as for such two-point closure models or the sub-grid-scale ones ( Sagaut & Germanol2005 ) which 
might either amplify or depress the thermalization effects. In the helical case, for instance, arti¬ 
ficial treatments (smoothing, hyperviscosity etc.) may strengthen the vortical mode and possibly 
unphysically (larger) conversion of acoustic energy to vortical energy, leading to underestimating 
the noise level. The techniques, such as the Kraichnan transformation and the helical representa¬ 
tion in the sharpened Helmholtz decomposition, in the treatment here may be applied to further 
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the extensions of|Leel(ll952h by Frisch et al. ( 1975 ). [Servidio. Matthaeus & Carbone (2008) and 
Zhu, Yang & Zhu[ ( 2014i to compressible magnetohydrodynamics, for similar discussions of the 


isotropic polarization issue etc. 

Concerning fundamental statistical physics, Eq. (I2.71 i provides some basis and motivates thoughts 
for the possibility of understanding in a more quantitatively definite way the polarization dis¬ 


cusse d here, by extending the well-known She-Jackson and Gallavotti-Cohen approaches (iGallavotti 


2014 and references therein). The idea is to replace the damping parameters in different compo¬ 


nents of the Navier-Stokes equations with some time-dependent ‘dynamical eddy-viscosity’ co¬ 
efficients (such that the system is also time-reversible, but, unlike the Galerkin-truncated invisid 
system, does not satisfy the Liouville theorem), to respect more observable(s), such as the helicity 
contained in subsets, say, each shell (She-Jackson), or even the whole set of fc (Gallavotti-Cohen), 
besides the energy. The Gallavotti-Cohen approach needs further theoretical assum ptions, such 
as the ‘chaotic hypothesis’ and ‘equivalent-ensemble conjecture’ (IGallavotti 112014t) : Such work 
for the simpler case of incompressible flow is formally easier to handle and is underway. Like 
the absolute equilibrium calculation, it may not be impossible that quantitatively different fluc¬ 
tuation relations for polarized components of the field could be established to definitely measure 
the compressible turbulence (isotropic) polarization. 
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